, Ground state energy of large atoms in a self-generated 
§ : magnetic field 

(N 

'—I 

, ■ Laszlo Erdos * 

Institute of Mathematics, University of Munich 
Theresienstr. 39, D-80333 Munich, Germany 

^ ^ and 

Jan Phihp Solovej 

fH ; Department of Mathematics, University of Copenhagen 

! Universitetsparken 5, DK-2100 Copenhagen, 

Denmark 



OO 



March 10, 2009 



Abstract 
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' We consider a large atom with nuclear charge Z described by non-relativistic quan- 

. turn mechanics with classical or quantized electromagnetic field. We prove that the 
absolute ground state energy, allowing for minimizing over all possible self-generated 

ij - electromagnetic fields, is given by the non-magnetic Thomas- Fermi theory to leading 

. order in the simultaneous Z — > oo, q — > limit if Za^ < k for some universal k, where 

' a is the fine structure constant. 
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1 Introduction 



The ground state energy of non-relativistic atoms and molecules with large nuclear charge 
Z can be described by Thomas- Fermi theory to leading order in the Z oo limit [Ll ILSj . 
Magnetic fields in this context were taken into account only as an external field, either a 
homogeneous one |LSYll[LSY2] or an inhomogeneous one [ESj but subject to certain regularity 
conditions. Self-generated magnetic fields, obtained from Maxwell's equation are not known 
to satisfy these conditions. In this paper we extend the validity of Thomas-Fermi theory by 
allowing a self-generated magnetic field that interacts with the electrons. This means, we 
look for the absolute ground state of the system, after minimizing for both the electron wave 
function and for the magnetic field and we show that the additional magnetic field does not 
change the leading order Thomas-Fermi energy. Apart from finite energy, no other assumption 
is assumed on the magnetic field. 

The nonrelativistic model of an atom in three spatial dimensions with nuclear charge Z > 1 
and with N electrons in a classical magnetic field is given by the Hamiltonian 




acting on the space of antisymmetric functions /\^ Ti. with a single particle Hilbert space 
Ti. = L^(R^) ® C^. The coordinates of the N electrons are denoted by x = (xi,X2, . . . ,xiy). 
The vector potential A : — > generates the magnetic field B = V x A and it can be 
chosen divergence-free, V ■ A = 0. The last term in (11.11) is the energy of the magnetic field. 
The kinetic energy of an electron is given by the Pauli operator 

T(A) = [o- ■ (p + A)]2 = (p + A)2 + ■ B, p = -zV. . 

Here cr is the vector of Pauli matrices. We use the convention that for any one-body operator 
T, the subscript in Tj indicates that the operator acts on the j-th variable, i.e. Tj{A) = 
[(Tj ■ {~iVx- + ^(xj))]^. The term —Z\xj\~^ describes the attraction of the j-th electron to 
the nucleus located at the origin and the term | the electrostatic repulsion between 

the i-th and j-th electron. 

Our units are f^{2me^)~^ for the length, 2me^Tr'^ for the energy and 2mech~^ for the 
magnetic vector potential, where m is the electron mass, e is the electron charge and h is 
the Planck constant. In these units, the only physical parameter that appears in (11.11) is 
the dimensionless fine structure constant a = e^(^c)~^. We will assume that Za"^ < n with 
some sufficiently small universal constant k < 1 and we will investigate the simultaneous limit 
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Z — > oo, a ^ 0. Note that the field energy is added to the total energy of the system and by 
the condition V • A = we have 

[ B^= [ |V®A|2, (1.2) 

where V® A denotes the 3x3 matrix of all derivatives diAj and | V® Ap = Yl^j=i I'^j^iP- We 
will always assume that the vector potential belongs to the space of divergence free if^-vector 
fields 

A:^{Ae H^-R^, R3), V • a = 0}. 

In the analogous nonrclativistic model of quantum electrodynamics, the electromagnetic 
vector potential is quantized. In the Coulomb gauge it is given by 

Aa(x) = A{x) = A_(x) + A_{x)* 

with 

Here g{k) is a cutoff function, satisfying \g{k)\ < 1 and supp g C {/c G : \k\ < A} with a 
constant A < oo (ultraviolet cutoff). The field operators A{x) depend on the cutoff function 
g{k) whose precise form is unimportant the only relevant parameter is A. For each k, the 
two polarization vectors e_(A;), e+(A;) G are chosen such that together with the direction 
of propagation k/\k\ they are orthonormal. The operators ax{k), ax{k)* are annihilation and 
creation operators acting on the bosonic Fock space over L^(R^) and satisfying the canonical 
commutation relations 

[ax{k),ay{k')] = [ax{k)\ay{k'r] = 0, [ax{k),ay{k'r] = Sxx'5{k - k'). 

The field energy is given by 

Hf = a-^ [ \k\y" ax{kyax{k)dk. 

The total Hamiltonian is 

<z - E [^.(Aa) - ut] + E u:^^ + ^/ (1-3) 

j=i I -'I i<j ' ' 

and it acts on (/\^ H) 
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The stability of atoms in a classical magnetic field [0 IFLLt ILLl ILLS] implies that the 
operator fll.ip is bounded from below uniformly in A. If ZoP' is small enough. It is known 
[LY[ IES2] that stability fails if Za^ is too large. The analogous stability result for quantized 
field |BFG] states that (11. 3p is bounded from below if Zo^ is small. In particular, we can 
for each fixed A define the operators in fll.ip and (11. 3p as the Friedrichs extensions of these 
operators defined on smooth functions with compact support. 

The ground state energy of the operator with a classical field is given by 



N 



E]^,^(A) = inf{(vI/,iJ^^^(A)M/) : G /\ (Co-(R^) ® C^) , ||^|| = l}, 

1 

and after minimizing in A we set 

E^,^= mf^E^,^(A). 

We note that it is sufficient to minimize over all A G where = H^(R.^,T{.^) denotes 
the space of compactly supported vector fields. It is easy to see that the Euler-Lagrange 
equations for the above minimizations in \& and A correspond to the stationary version of 
the coupled Maxwell-Pauli system, i.e., the eigenvalue problem if^^(A)\E' = E'^ together 
with the Maxwell equation V x B = Ana'^J^f,, where J<i, is the current of the wave function 
It is for this reason that it is natural to refer to B as a self-generated magnetic field in this 
context. 

In the case of the quantized field, we define 

N 

e^'^^ = m[{^,h'^';'^^) : vi/g/\(Co^(r3)®c2)®^, ||^|| = i}. 

1 

The stability results of [B iFLLl iLLl ILLSI IBFGJ imply that > -oo and E^^^^ > if 
Za^ is small enough. 

Finally, we define the ground state energy with no magnetic field as 



En,z ■= = 0). 



In all three cases we define 



Et := inf^E*^^, # e {cl, qed, nf}, 

for the absolute (grand canonical) ground state energy. The main result of this paper states 
that the magnetic field does not change the leading term of the absolute ground state energy 
of a large atom in the Z —>■ (yo limit. In particular, Thomas- Fermi theory is correct to leading 
order even with including self-generated magnetic field. 
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Theorem 1.1. There exists a positive constant k such that ij Zo? < k, then 



(1.4) 



For the quantized case, if we additionally assume 



1 _7 _ 1 



(1.5) 



with some < 7 < then 



+ CZak^ > Ef^ > Ef - Czi-\ 



(1.6) 



We note that Zak^ < Z'^''^ «/A < ^-1/4^11/12 the Z ^ 00 limit. 

Remark 1. The leading term asymptotics of the non- magnetic problem is given by the 
Thomas- Fermi theory and E"^ = — ctf^^ + O(Z^) as Z ^ 00, where ctf = 3.678 ■ (37r^)3 
is the Thomas- Fermi constant. The leading order asymptotics was established in [LSj (see 
also [L]). The correction to order is known as the Scott correction and was established in 
[Hi ISWlj and for molecules in [IS] (See also [SW2l[SW3llSSj ). The next term in the expansion 
of order Z^^^ was rigorously established for atoms in [FSj . 

Remark 2. The exponents in the error terms are far from being optimal. They can be 
improved by strengthening our general semiclassical result Lemma 11.31 for special Coulomb 
like potentials using multi-scale analysis. 

Remark 3. For simplicity, we state and prove our results for atoms, but the same proofs 
work for molecules as well, if the number of nuclei K is fixed, each has a charge Z, and assume 
that the nuclei centers . . . ,-Re'} are at least at distance away, i.e. \Ri — Rj\ > 



Remark 4- Theorem 11.11 holds for the magnetic Schrodinger operator as well, i.e. if we 
replace the Pauli operator T(A) = [cr • (p + A)]^ by T(A) = (p + A)^ everywhere. The proof 
is a trivial modification of the Pauli case. The argument is in fact even easier; instead of 
the magnetic Lieb-Thirring inequality for the Pauli operator one uses the usual Lieb-Thirring 
inequality that holds for the magnetic Schrodinger operator uniformly in the magnetic field. 
We leave the details to the reader. Note that although the condition Za"^ < k is not needed in 
order to ensure stability in the Schrodinger case, we still need it in the statement in Theorem 
11.11 In the Schrodinger case this condition is not optimal. 

The upper bound in (11.41) is trivial by using a non-magnetic trial state. The upper bound 
in (II. 6p is obtained by a trial state that is the tensor product of a non-magnetic electronic 
trial function with the vacuum \Q) of JF. The field energy Hf and all terms that are linear in 



cZ-y^i^j. 
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A give zero expectation value in the vacuum. The only effect of the quantized field is in the 
nonlinear term A^. A simple calculation shows that 

{n\A^\Q) < Cak^. 

The main task is to prove the lower bounds. Using the results from [BFG] . the result for 
the quantized field (11.61) will directly follow from an analogous result for a slightly modified 
Hamiltonian with a classical field. Let 

N 

i^^,z(A)=i7^,^,„(A) = 5:[T.(A)-A]+J^^^ + J_ / iv^Ar (1.7) 
with some 

r = DZ-^'^ with D>1. (1.8) 

Note that instead of the local field energy, the total local ^-norm of A is added in (11.71) . By 
(II. 2p . we have 

Hf,^,{A) > H^^A) (1.9) 
for any A ^ A. We define the ground state energy of the modified Hamiltonian (11.70 

AT 

i?^,^(A) :=inf |(vI/,i7^,^(A)vl/) : ^i! e /\ {C^ {K') C) , ||vl>|| = l}, 

1 

and set 

En,z ■= inf En,z{A), Ez := mfEN,z, 

where the infimum for A G ^ can again be restricted to compactly supported vector potentials 
A E Aq. For the modified classical Hamiltonian we have the following theorem: 

Theorem 1.2. Let Zo? < k and assume that r = DZ'^/^ with 1<D< Z^/^'\ Then 

Ez >Ez> Ef - CZ^/^D-\ (1.10) 

Taking into account (II. 9p , Theorem 11.21 immediately implies the lower bound in (II. 4p . The 
proof of the lower bound in (11.60 follows from Theorem 11.21 adapting an argument in |BFGj 
that we will review in Section [6] for completeness. 

One of the key ingredients of the proof of Theorem 11.21 is the following semiclassical state- 
ment that is of interest in itself. The first version is formulated under general conditions but 
without an effective error term. In our proof we actually use the second version that has a 
quantitative error term. 
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Theorem 1.3. Let Th{A) = [cr ■ {hp + A)]^ or Th{A) = {hp + Af, h<l, andV > 0. 

1) IfV e L5/2(r3) nL4(R3), then 

TT[Th{A)-V]_ + h-'^ [ >Ti[-h^A-V]_+o{h-^) as h 0. (1.11) 

2) Assume that ||V^||oo ^ K with some 1 < K < Ch^"^ and consider the operators with 
Dirichlet boundary condition on Q G R'^. Let Bji denote the ball of radius R about the origin 
and let := ^2 + denote the y/h-neighborhood of the set ^l. We set \^^\ for the 
Lebesgue measure of fl^. Then 

Tr [(T,(A) - V)n] + h'^ [ 

Jn^ (1.12) 

> Tr [{-h'A - V)n]_ - Ch-'K''/'\Q^\{hK'/^Y^' [l + [hK'/^Y^']. 

Remark. Despite that the electrons are confined to Q, their motion generates a magnetic field 
in the whole R^, so the magnetic field energy in (11.121) is given by integration over R^. 

We use the convention that letters C, c denote positive universal constants whose values 
may change from line to line. 

2 Reduction to the Main Lemmas 

Proof of Theorem li.M We focus on the lower bound, the upper bound is trivial. We start 
with two localizations, one on scale r > and the other one on scale d < Z'^/^. The 

first one is designed to address the difficulty that the if^-norm of A is available only locally 
around the nucleus. This step would not be needed for the direct proof of ( 11.40 . The second 
localization removes the "Coulomb tooth", i.e. the Coulomb singularity near the nucleus. 

In this section we reduce the proof of the lower bound in Theorem 11.21 to two lemmas. 
Lemma [2TT] will show that the Coulomb tooth is indeed negligible. Lemma [2.21 shows that the 
magnetic field cannot substantially lower the energy for the problem without the Coulomb 
tooth. In the proof of Lemma 12.21 we will use Theorem 11.31 

Recall that denotes the ball of radius R about the origin. We construct a pair of 
smooth cutoff functions satisfying the following conditions 

el + ef = 1, supp 01 C B2d, Oi = l on Bd, iV^ol, |V^i| < Cd'K 

We will choose 

d = 6Z-^/^ (2.13) 
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with some S < 1, in particular d < r. 
We split the Hamiltonian as 

with 



N 



i=l 



eo(T{A)- — -{\veo\' + \ve,\'))eo 
1 1 



i<j 



N 



1=1 



h(T{A) 



Z 

\x\ 



IV® A| 



(2.14) 



{\veo\' + \vei\'))ei 



+ 



Bar 



where we used the IMS localization formula that is valid for the Pauli operator as well as for 
the Schrodinger operator. 

In Section [3] we deal with Hj^ ^, to prove that it is negligible: 

Lemma 2.1. There is a positive universal constant k, such that for any Z,a with Zo? < k 
we have 

inf inf H],z{A-) > -CZ'/'^S'^' - Z^^'S-' 

N AeAo 

if CZ^'^^^ < 5 < D with a sufficiently large constant C. 

Starting Section H] we will treat H'^ ^{A) and we prove the following 

Lemma 2.2. There is a positive universal constant n such that for any Z, a with Zo? < k 
we have 

m{MH%z{-A) > -ctf^^/^ - CZ^/3[Z-i/30 + 2^-1] (2.15) 
with a sufficiently large constant C if Z^^l'° < 5 < 1 and D < Z^l'^^b^^l^'° . 

The main ingredient in the proof is Theorem 11.31 that will be proven in Section O The proof 
of the lower bound in Theorem 11.21 then follows from Lemmas 12.11 and 12.21 after choosing 

5 = Z-2/63. □ 
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3 Estimating the Coulomb Tooth 



Proof of the Lemma \2.1\ Let xo be a smooth cutoff function supported on B^r sucli tliat 
|Vxo| < Cr^'^ and xo = 1 on B2r- Let (A) := l-Bsr]"^ /g^ A. We define 



Ao := (A - (A))xo, Bq := V x Aq, 
tlien V (g) Ao = XoV (g) A + (A - (A)) (g) Vxo- Clearly 

/ / |V®AoP<2/ Xo|V® Ap + Cr-2 / (A - (A)) 



(3.16) 



< Ci / |V® A|2 



(3.17) 



for some universal constant Ci, where in the last step we used the Poincare inequality. Let 
be a real phase such that V(f = (A). Since Xo = 1 on the support of Oi by D > 6, we have 

0iT(A)0i = eie-'^T{A - {A))e'^ei = 0ie-^^T(Ao)e'^^i. 

After these localizations, we have 

N 



i=i L 



^ie-*'^(T(Ao) -iy(x))e*% 



1 



2 / ^0 



(3.18) 



with 



W{x) 



cd- 



\x\ 



l{\x\ < 2d) . 



Now we use the "running energy scale" argument in |LLSj . 



N 

E 



T(Ao) - W 



> 



M.e{T{Ao) - W)de 



n jj^ /*00 

>-J U-e{T{Ao)-W)de- j Aro(^T(Ao) -ly + e)de (3-19) 
>- /V„e(T(Ao)-H^)de- rAro(T(Ao)--l^ + -)de, 



where M-e{A) denotes the number of eigenvalues of a self-adjoint operator A below — e. 
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In the first term we use the bound T{Aq) > (p + Aq)^ — |Bo| and the CLR bound: 



r U-e{T{AQ) - W)de <C r de f {W + |Bo| - e] 

Jo Jo JR3 



3/2 



<C I del {W -e/2fl'^ + C I del (|Bo| - eV2^) 



R3 



\3/2 



R3 



JR^ Jrs 

= ^ ^^-2 _^ ^^1/2 /■ 



(3.20) 



In the second term of f l3.19p we use 
1 



T(Ao) 
and that 



-W > 
/i - 2 



(P + Ao)' 



2eZ 



1(1x1 < 2d) > (p + Ao 



[P + Ao) 



iBr 



Ce 



l{\x\ < 2d), 



4eZ /2eZ\2 



> 



(3.21) 



I.e. 



T(Ao)--W^> -(p + Ao)'-2 

jj, 2 \ jj, 



} -|Bc 



Ce 



l(|x| < 2d). 



We choose /x = 4Z^, then using Ce/ ^.d'^ < e^/4/i for /i < e (i.e. C < ((5Z^/^)^), we get 



AAo(T(Ao)--iy + -)de< / Aro(-(p + Ao) 



IBr 



4/x 



de 



<C Tde [ (|Bo|-eV4/x) 

Jo ^R3 

<C/.V^/ Bg. 

JR3 



(3.22) 



Note that if Za"^ < k, with some sufficiently small universal constant k, then the magnetic 
energy terms in (13.201) and (I3.22p can be controlled by the corresponding term in (I3.18p . 
Combining the estimates (KIM . (KW) . (ICTD and (Km we obtain 



Hhzi-^) > -CZ^/^d^l^ - Cd- 



(3.23) 



and Lemma [2.11 follows. □ 
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4 Removing the magnetic field 



Proof of Lemma \2.^ We start with two preparations. In Section 14.11 we give an upper bound 
for the number of electrons in the truncated model described by ^{A). In Section 1^2] 
we then reduce the problem to a one-body semiclassical statement on boxes. The semiclassical 
problem will be investigated in Section [5] and this will complete the proof of Lemma 12.21 



4.1 Upper bound on the number of electrons 

Let 

N 

1 

be the ground state energy of the truncated Hamiltonian if^^(A) defined in (12.141) . The 
following lemma shows that we can assume N < CZ when taking the infimum over N in 
fl2.15p . The proof is a slight modification of the proof of the Ruskai-Sigal theorem as presented 
in |CFKSj . We note that the original proof was given for the non- magnetic case and it can 
be trivially extended to the Schrodinger operator with a magnetic field but not to the Pauli 
operator. This is because a key element of the proof, the standard lower bound on the hydrogen 
atom, —A — Z/\x\ > —Z'^/A, is valid if —A = replaced by (p + A)^ but there is no lower 
bound for the ground state energy of the hydrogen atom with the Pauli kinetic energy that is 
independent of the magnetic field. However, for the truncated Coulomb potential the trivial 
lower bound can be used. 

Lemma 4.1. There exist universal constants c and C such that for any fixed A ^ Aq and Z 
we have 

whenever N > CZ and Z^^l'° < 5 < c. In 'particular 



inf inf E^7(A) = inf lYvi E% yiA) (4.24) 



%fZ-^l^ <5<c. 



Proof. We mostly follow the proof of Theorem 3.15 in |CFKSj and we will indicate only 
the necessary changes. For any x = {xi,X2, ■ ■ - Xn) G R'^^ we define 

Xoo{x) := max{|xj|, : i = 1,2, . . . , N} 

Ao := {x : | < ^ V j = 1, 2, . . . A^} 

Ai := {x : \xi\ > (1 — C)xoo{x), Xoo{x) > — } 
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for some fixed positive q and C, < 1/2 to be chosen later. According to Lemma 3.16 in |CFKSj . 
there is partition of unity { Ji}i=o,i,...Ar, with J"^ = 1, supp Ji C Ai such that the gradient 
estimates 

L{x) = J2 |VJ^(a;)P < — if x e 



i=0 
N 



L{x) = |VJ.(x)|2 < — if X G Aj, J > 1 



1=0 



hold with a suitable universal constant C. Moreover, Jq is symmetric in all variables, while 
Jj, i > 1, is symmetric in all variables except Xj. 

We subtract the local field energy that is an irrelevant constant, i.e. define 

and Em = inf Spec H^. We will show that En = E^^i for > CZ. By removing one 
electron to infinity, clearly Ej^ < E^^^i < 0; here we used the fact that A is compactly 
supported. 

By the IMS localization 

N 

Hn = MHn - L)Jo + UHn - L)J,. (4.25) 

i=l 

In the first term we use that on the support of we have —Z\x\~^ > —Zd^^. Hence 

Jo{Hn - L)Jo >Jo(- CZNd-' - CNd~^ + —— ^ ^) Jq- (4.26) 

Choosing ^ = 8rf we see that Jo{Hn - L)Jo> if N > CZ with a constant C if 6 > CZ-'^I'^. 
To estimate the terms Ji^Hn — L)Ji for i 7^ 0, we define 



TV 



9„(r(A) - ^ - (iv9„r^+ iv9,p))9„ + y] ^ 

V \x\ / ^-^ \Xk — Xi\ 
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On the support of Jj we have \xi\ > g/A = 2d, so V^o and V^i vanish. Then we can estimate 



Ji{HN — L)Ji >Ji I i^^Li — -j — r + 



- 1 CN^/^ 



Ji 



>Ji I En-i + -j — j- 



^^i^-O-z 



J,; 



(4.27) 



a N>CZ and is large. Thus we conclude from Km . Km and Km that En > En-i 
if > CZ. D 



4.2 Reduction to a one-body problem 

We start by presenting an abstract lemma whose proof is given in Appendix Rl 

Lemma 4.2. Let ^ he a one-particle operator on Ti = L^(R^) and let W he a two-particle 
operator defined on HAH. We assume that the domains o/f) and W include the functions. 
Let 9 e C°°(R'^) with compact support VL := suppO. Then 



N 



N 



inf 



l<i<j<N 



N\ 



1^1 



i=l 



(n n 
(•^'(^^5^+ E ^■^)^) ■ '^e/\Co-(^]), ||$|| = l 
i=l l<i<j<n 1 



(4.28) 



where f)j denotes the operator [) acting on the i-th component of the tensor product, and similar 
convention is used for the two-particle operators. The same result holds with ohvious changes 
ifH = L\R^)®C\ 

To continue the proof of Lemma l2.2| we first localize if^^(A) onto a ball Bj. of radius 
r = DZ^^^^ (see ( 11.81) ) and we also localize the magnetic field as in Section [3J We introduce 
smooth cutoff functions Xo and Xi with 



xl + xl = 1, supp xo C B2r, Xo = ^ on B,., | Vxol, |Vxi| < Cr ^ 
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We get 



N 



oXoe 



1=1 



+ 



r,(Ao) - -- 

Xi . 



i<j 



(4.29) 



R3 



using that the new locahzation error |VxiP + |VxoP ^ Cr~^ < Cc?"^ and that —Z/\x\ > 
—Zr~^ on the support of Xi- We also used fl3.17p . 

Let Ad^r = {x : d < \x\ < r} G R^. Using (14.241) . the positivity of the Coulomb repulsion 
\xi — Xj\~^ > and Lemma 112] with 9 := ^^oXo we obtain 



inf inf H%JA) > inf 

N Aey^o ' N<CZKQ(iAo 



\ 1=1 l<] 



+ 



R3 



-1 



(4.30) 



where the infimum is over all antisymmetric wave functions ^ e C^{Ad,r) ® with 
||\1/||2 = 1. The notation [H]q indicates the A^-particle operator H with Dirichlet boundary 
condition on the domain C R'^^. 



We define 



f{x)g{y) 

R3xR3 — y\ 



dxdy. 



Lemma 4.3. There is a universal constant Co > such that for any 

with II II 2 = Ij /or any nonnegative function g : R^ R wi/i -D(^, g) < oo, for any A G 

and for any e > we have 



N 



vI/,U5^T.(A)+5^ 



i=l 



X i j I " 



Cn 



R3 



N 



(4.31) 



> -D{g, g) + [^,J2{Q 



* \x. 



1-1 



1=1 
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N , (4-32) 

4/3 



Proof. By the Lieb-Oxford inequality [LOj and by the positivity of the quadratic form 

N 

> -D{g, g) + U, y^ig * \xr')^) - C / g'^ 

where gis,{x) is the one-particle density of \E'. 

The error term is controlled by the following kinetic energy inequality for the Pauli operator 

N 

(v[/, [gT(A),]vI/) >cy^^ min{^,f, 7//} -7 (4.33) 

with some positive universal constant c and for any 7 > 0. For the proof of (14.331) use the 
magnetic Lieb-Thirring inequality 

N 

[V[T(A) -f/]J^\ > -C /" U'^'^-C^-^l U^-if B^. 



With the choice U = (3mm{g^^ ,'yg]^^} we can ensure that ^Ug^, > CU^^"^ + C'j'^U'^ if f3 i 
sufficiently small (independent of 7) and this proves (14.331) . 
Thus 



IS 



R3 



4/3 ^ 



"M min{^^/^7^^''^} + 7 / g^ 

JR3 JR3 

N „ 

< (c7)-i(^, [5^r(A),]^)+c-iy^^B2 + 7iV 



(4.34) 



so choosing 7 = Ce ^ with a sufficiently large constant C, we obtain (I4.3ip . □. 

Using Lemma 14.31 we can continue the estimate (14.301) (with writing A instead of Aq in 
the infimum) as 

infinfi7^^(A)>(l-£) inf inf <^ inf [ V[T(A) + W^] J A + 

- D{g, g) - Ce-^Z - CZ^'^d''^ - CZ'^'^D'^ 

(4.35) 
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with 



W{x) :-- 



1-e 



— -r-r + Q * \x 

\x\ 



-1 



and assuming that a < with some small universal Oq. 

We now perform a rescaling: x = Z^^^'^X, p = Z^^^P and 

A(X) = Z-2/3a(Z-i/=^X), B(X) = V X A{X) = Z~^B{Z-^/^X) 
Introducing h = Z"^!"^ and T/i(A) := [(/iP + A)-cr]^, we obtain that the kinetic energy changes 



as 



[(p + A) ■ cr]2 = Z^!%Z-^^''V + A) ■ erf = Z^/^Th{A.) 



and the field energy changes as 



B2(x)da; = Z / B^{X)dX. 



R3 



R3 



The new potential energy is 



w{x) = z-^i'^w{z-^i^x) 



where q{X^ = Z '^g{Z ^/^X) and D{g, g) = Z '^^^D{g, g). After rescaling, we get from (14.351) 



inf inf Hi J A) 



N 



> (1 - inf inf <^ inf V [T;,(A) + W], 

- Z'/^'Dig, g) - Ce-^Z - CZ^''^-^ - CZ'/'D-\ 



\& ) H 

As,D I CZa^ J^3 



(4.36) 



where As^d = {X : S < \X\ < D} and inf,j, denotes infimum over all normalized antisym- 
metric functions. Using ( 11.12^ from Theorem 11.31 and the fact that Za"^ < n, we get 



inf inf i70 ^^(A) > (1 - e)Z^^''TT [{-h^A + W)asJ- - ^^^3/6^3^-13/4 
- Z''/''D{g, g) - Ce-'Z - CZ^'IH-'' - CZ'^'^D-^ 



(4.37) 



assuming 5 > Z By standard semiclassical result for Coulomb-like potentials (see e.g. 
the result in Section V.2 of jL]): 

Tr [(-/I'A + W)a,^]- > Tr [ - /i^A + W]^ > -C^ch-^ [ Wl^^ - Ch-^+^'^l (4.38) 



16 



where Csc = 2/(157r^) is the Weyl constant in semiclassics. The ^ exponent in the error term 
is far from being optimal; the methods developed to prove the Scott correction can yield an 
exponent up to one (see Remark 1 after Theorem ll.il) . 

Taking the optimal ^to be the Thomas- Fermi density for Z = 1 ^ = ^^tf (see, e.g. Section 
II of [L]) and defining the Thomas- Fermi constant as 



Ctf := D{gTF, Qtf) + C, 



sc 



R3 L 



5/2 



we get 



inf inf H% y 



( 



1 

Jxi 



g* \X 



-1 



5/2 



- CZ^/3-1/30 _ ^^13/6^3^-13/4 _ ^^-1^ _ CZ5/3^-2 _ CZ'^^D-' 
> _(1 _ 5)-3/2cTp^7/3 _ ^^^7/3-1/30 _ ^^-1^ _ CZ'/'^d'^ - CZ'/'^D^' 

(4.39) 

where we optimized for e and we used that D < Z^^'^^5^^^^^ and Z~^^^ ^ S < 1. This completes 
the proof of Lemma 12.21 □ 



5 Semiclassics: proof of Theorem 11.3 



We present the Schrodinger and Pauli cases in parallel. We prove the statement with Dirichlet 
boundary conditions fll.l2p in details and in Section [5^ we comment on the necessary changes 
for the proof of the (11.111) . The potential V is defined only on Q, but we extend it to be zero 
on R'^ \ we continue to denote by V its extension. 

5.1 Localization onto boxes 

We choose a length L with h < L < ^h^^"^- Let = + Bi be the L-neighborhood of f2. Let 
Qj, = {y eR^ : \\y - k\\oc < L/2} with k E {LZf n denote a non-overlapping covering 
of with boxes of size L. In this section the index k will always run over the set {LZy^nni. 
Let be a partition of unity, J2k^k = 1' subordinated to the collection of boxes Qk, such 
that 

supp C (2Q)fe, iV^fcl < 
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where {2Q)k denotes the cube of side- length 2L with center k. Let be a cutoff function such 
that = 1 on {2Q)k (i.e. on the support of ^fc), supp^^ C Qk '■= (3(5)^ and |V^fe| < CL~^. 
Let (A)fc = IQfcl""'^ In Similarly to fl3.16p . we define := (A — {A)k)^k and Bjt : = 



V X At., then 



Bl<C |V® A| 

R3 Jq^ 



(5.40) 



as in fl3.17p . From the IMS localization with ipk satisfying hVipk = (A^) we have 



Tr 



[ThiA) - V]n\ + /" = inf Tr (i[Th{A) - V]) + h'^ [ | V ® A 



>inf Yl ^^(^) 



(5.41) 



fce(LZ)3nQi 



with 



Sk{7) := Tr 7ae-*'^'=[T;,(A - (A),) - Vy^'^k - l\hV^k 



+ Coh~ 



IV® A| 



with some universal constant cq. Here inf* denotes infimum over all density matrices < 
7 < 1 that are supported on fl, i.e. they are operators on L^(f2) ® C^. We also used 
/j^3 = Jj^3 I V (S> Ap and we reallocated the second integral. We introduce the notation 



J^k ■= Coh~ 



IV® Ar. 



5.2 Apriori bound on the local field energy 

In case of the Pauli operator, for each fixed box Qk we apply the magnetic Lieb-Thirring 
inequality [LLSj together with (I5.40p and with the bound ||V^||oo < K to obtain that for any 
density matrix 7 



^fc(7)>Tr [mAk)-V-Ch'L 



2 r-2i 



> -Ch-' [ [V + Ch''L-^f^~c( [ [V + Ch''L-^]^y(h-^ [ B^)'^\J'^ 



> -c 



2 



V® A M + J-fc 



Qk 



(5.42) 
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using h < L and 1 < K < Ch~'^. In the Schrodinger case we use the usual Lieb-Thirring 
inequahty jLTj that holds with a magnetic field as well. The estimate fl5.42l) is then valid even 
without the third term in the second line. 

Let 5* C {LTi)'^ n fix denote the set of those k indices such that 



holds with some large constant C. In particular 

^k{l) > 0, for all A; ^ S" and for any 7. 

5.3 Improved bound 

We use the Schwarz inequality in the form 

niA - (A)fc) > -(1 - ek)h''A - Ce,\A - {A),y, 
with some < < |- We have for any 7 supported on Q that 

Skii) >Tr - 2ek)h'A -V - Ch''L-^]ikh 



+ Tr 



l.\-euh'A-Ce,\A-{A)uY]l. 



We will show at the end of the section that 



Tr 



> Tr 

Using flCTj) and ^M), 

k keS 

> 5^ Tr [inU-h^A - V]^kln] _ + Y.^^ 

k kes 

> inf Tr [^7^61^1-/^' A - V]ln \ 

k 

>Tr [(-/i2A-V)f,]_ + ^P, 



(5.43) 



(5.44) 



(5.45) 



(5.46) 



(5.47) 



kes 
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with 



Vk := Tr 



In the last step in (15.471) we used that for any collection of density matrices 7^, the density 
matrix IcCfcTfc^fclf^ is admissible in the variational principle 

Tr [(-/i^A-r)n]_ = inf |Tr7[-/i2A- V] : < 7 < 1, supp 7 C ^]}. (5.49) 

We estimate T>k for A; G as follows 

Vk>-Cel^h~'' [ iA-{A)k)'-Ch-'K'^/^\Qk\{ek + h^L-^)+J'k 

Jq, (5.50) 

In the first step we used Lieb-Thirring inequality, in the second step Holder and Sobolev 
inequalities in the form 



/ (A - (A),)^ < CL'/' ( / I V ® Ap) 



5/2 



We choose 

and using the apriori bound (I5.43p . we see that 

Sk < Ch-'/^LK'/\ 

Thus, assuming 

L < ch'/^K~^/^ (5.51) 

with a sufficiently small constant c, we get Ek < 1/3. With this choice of Ek, and recalling 
iQfcl = Q\Qk\ = 9L^, we have 

If we choose L = h^/'^K~^/^, then 

1/2 
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This choice is allowed by fl5.5ip li K < ch li ch < K < h then we choose 
L = ch^/'^K"^/'^ and we get from fl5.52p 



Combining these two inequalities, we get that 

Vk > -Ch~^L'K'>'\hK^/^Y'\l + [hK^/^y^] 
always holds. Summing up fl5.53p for all k and using that 

L'<C\^,l\<C\^^\ 



(5.53) 



(recall that is a neighborhood of and 3L < h^^'^), we obtain from fl5.47p and ([51 



irif 5^^,(7) >Tr [{-h'A - V)n\_- Ch-'K'/'\n^,\{hK'/y'[l + {hK^/y'] (5.54) 

k 

and this proves f ll.l2p . 



Finally, we prove (I5.46p . Let 7 be a trial density matrix for the left hand side of (I5.46p . 
We can assume that 



> Tr 



Then 
>Tr 

>Tr 



1 







Tr 



llnU-7:h'^ -V- Ch'L~' - K]^^^ 




lln^k[--h''A + K]iuln - Ch-' [V + K + Ch^L-^f\ 
6 J Jq^ 



(5.55) 



where we used Lieb-Thirring inequality. Thus, using \V\ < K , h < L and K > 1, we have 



Tr 



o J 



Therefore 



Tr 



7lnefc[-(l - 2ek)h'A -V- Ch'L-']^kl 



> Tr 



(5.56) 



Now f l5.46p follows by variational principle. □ 
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5.4 Reduction of ( TrTTj) to ( TTT^ 



We approximate V G L^/^ fl L'^ by a bounded potential ||V^||oo < that is supported on 
a ball -Bi?/2 and V <V. By choosing and R sufficiently large, we can make \\V — V\\5/2 + 
\\V — V\\4 arbitrarily small. We choose a cutoff function xr that is supported on Br, xr = ^ 
on Bji/2 and |Vxi?| < CR^^ and let xr satisfy x'r + Xr = 1- 

Borrowing a small part of the kinetic energy, by IMS localization we have 

n{A) - V >a - e)xR[niA) - V]XR ^ 

+ eUA) - (V - (1 - e)V) - \Vxr\' ~ \Vxr\' 
Using the magnetic Lieb-Thirring inequality [LLSj to estimate the second term, we get 
Tr [T,(A) - >(1 - e)Tt [(T,(A) - V)n,]^ 

I,-, r ^2 (5-58) 



Ce-'/^h-' [ \uf^-cj \U\^--h-' [ 



with 



U:={V-{1- e)V) + \Vxr\' + \Vxr\' . 
For the first term in (15.581) we use (11.121) (and that it holds even with a 1/2 in front of /B^) 



and the fact that 



Tr [{-h^A - V)b^] _>TT[-h^A-V]_ 



by monotonicity, V < V. The second and the third terms in (I5.58P can be made arbitrarily 
small compared with for any fixed e if R and K are sufficiently large and h is small. 
Finally, choosing e sufficiently small, we proved (11.111) . □ 

6 Proof of the quantized field case 

For the proof of the lower bound in (II. 6p . we follow the argument of [BFG] to reduce the 
problem to the classical bound (11.101) . We set 



Hg = a-' [ \g{k)\^\k\y^ax{kyax{k)dk. 



to be the cutoff field energy, then Hf > Hg and only the modes appearing in Hg interact with 
the electron. By Lemma 3 of [BFGj . for any real function / G L^(R^) fl L°°(R^) we have 



^ I f{x)\V^A{x)\Mx<a^f\\^Hg + CaA\iJii,. 
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Applying it with / being the characteristic function of the ball B^^ with r = DZ (the 
radius of the ball is here chosen differently from |BFGj ) and using Za"^ < k we get 

Hf>Hg> (^)Hg |V ® A{x)\Mx - CK-^aM^D^. 

Setting a = (k/ZY^"^, i.e. Za"^ = k, we have for k sufficiently small 

where Ei\i^z,a is the ground state energy of the Hamiltonian (11.71) with fine structure constant 
a. Applying (ll.lOp to this Hamiltonian, we get 

Ef'' > Ef - CzId-^ - CK-^aA^D^ 

whenever 1 < D < Z^ . Writing D = Z"' and applying the upper bound (ll.Sp on A, we obtain 
the lower bound in (11.61) . □ 



A Proof of Lemma 4.2 



Let the function xi^) ^ C°^(R^) be defined such that 9'^{x) + x^i^) = 1- any subset 
a C {1, 2, ... , A^} we denote by Xa the collection of variables {xj : i G a} and define 



a • — •a\'^a 



:= n^(^^)- 



We set the notation a'^ = {1, 2, ... , A^} \ a for the complement of the set a and set n :- 
{1, 2, . . . , n}. Let |a| denote the cardinality of the set a. 

For an arbitrary function ^ G /\i C^(R^), ||\E'|| = 1, and for < < A^ we define 



r" := eJTr„c 



S^c|^')(^|S„c 



where Tr„c denotes taking partial trace for the Xn+i,Xn+2, ■ ■ ■ ,xn variables. Define the 



fermionic Fock space as JF = 0^^q7-^„ with 7i„ := /\^H and we define a density matrix 



r:= Yl rH=f;('^y" on ^. 
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We first prove that F < / on It is sufficient to show that F < 7 on the n-particle sectors 
for each n. Let n < N, choose $ e Tin, and compute 

aC{l,2,...iV} aC{l,2,...iV} ^ 

|a|=n |a|=n 

= / da;adx'„dya<=^(a;a)©a(2;a)Sac(yac)*(x«, yac)'l>{x'^, yac)E^c{y^c)e^{x'J^{x'J 

<J2 I dXadx'^dyacel{Xa)El.iyac)\^{Xa,yac)\''mx'J\'' 

-ml 

(A.59) 

using Schwarz inequahty and that 1 = njLi[^^(^j) +X^(^j)] — X^a ®Q;(^a)"Q;=(^Q;=)- 
Second, for a fixed n < N, we compute 

n=0 i=l n=Q ^ ^ i=l 

— y ] / dXadXac^iXa, Xac)Qa{Xa)^a''{Xa'')^i(^^a{Xa)^a''{Xa^)^{Xa, Xac)^ 

a iea 

(A.60) 

= ^ ^ / dyi'if{Xa,Xao)Ql\^q{Xa\{i})'E!^a-i^a-)d{Xi)i)i(9{Xi^^ 
1=1 a : iEa 

N 



1=1 



where the trace on the left hand side is computed on J^. In the last step we used that for 
any fixed i, we have 1 = Ilj^J^^C^j) + X^i^j)] — Z]a^l"l= where the summation is over all 
a C {1, 2, . . . , A^} \ {i} and S'^ = {1, 2, . . . , iV} \ {i} \ a. 

A similar calculation for the two-body potential shows that 

N 

n=0 l<i<j<n l<i<j<N 
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Thus, by the variational principle, 

N 



i=l i<i<j<N 



N n 



n=0 1=1 l<«<i<" 



Since F is a density matrix supported on fi, we obtain (14.281) . □ 
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